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_ _ /We  have  also  in  this  study  demonstrated  the  possibility  of  utilizing 

the  image  subtraction  technique  as  applied  to  micro-circuit  board 
inspection,  which  would  have  a  profound  effect  on  precision  inspection 
scheme  for  large  scale  automatic  assembly  for  U.S.  Army  needs  and  her 
supporting  industries.  This  technique  would  provide  the  capability  of 
rapid  identification  inspection,  and  possibly  utilizing  for  synthesis  and 
fabrication.  The  effects  due  to  coherence  requirement  are  also  included. 
Experimental  simulations  of  the  IC  chip  inspection  are  provided.  _  - 

We  have  also  reported  an  application  of  this  source  encoding  technique 
to  generate  a  2-D  color  hologram  with  an  extended  incoherent  source.  The 
color  hologram  image  can  be  retrieved  with  a  simple  white-light  optical 
processor.  Experimental  demonstration  of  a  color  hologram  image  generated 
by  this  technique  is  included.  We  have  noted  that  a  3-D  hologram  image 
for  phase-type  object  of  definite  depth  may  be  generated  by  this  technique 

Finally,  the  net  effect  is  to  emphasize  the  truth  of  our  initial 
assumption  of  the  white-light  image  subtraction  can  be  an  essential  part 
of  the  missile  tracking  and  identification  application.  This  research 
effort  in  real-time  tracking  program  should  be  continuing  to  pursue. 
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I .  Introduction 
1 . 1  U . S .  Army  Goals  and  Needs 

The  U.S.  Army  needs  to  improve  the  tracking  technology  not 
only  by  improving  the  simplicity  of  the  tracking  equipment  but 
also  needs  to  improve  the  quality  and  versatility  of  the  tech¬ 
nique. 

The  basic  tracking  technique  that  we  would  like  to  investi¬ 
gate  may  center  at  white-light  optical  image  subtraction  tech¬ 
nique;  particularly  in  the  application  to  missile  tracking  and 
identification.  We  would  like  to  stress  that  the  image  sub¬ 
traction  techniaue,  as  compared  with  the  optical  correlation 
detection  technique'*',  in  the  application  to  missile  tracking 
offers  the  following  advantages: 

1.  The  tracking  technique  is  generally  free  from  scale 
and  object  orientation  problems. 

2.  The  system  is  not  limited  to  a  small  number  of 
objects  to  be  detected. 

3.  Since  the  system  does  not  emplov  matched  spatial 
filter,  the  tracking  system  is  versatile  and  simpler 
to  operate. 

U.  The  system  is  generally  smaller  and  economical  to 
install . 

There  are  also  several  reasons  for  selecting  the  optical 
image  subtraction  technique  rather  than  their  electronic  or 
digital  counterpart,  because  the  optical  technique  has  the 
following  capabilities:  1 

1.  It  is  capable  of  handling  very  large  information 


elements . 


2.  It  has  a  very  large  space-bandwidth  product. 

3.  It  has  the  capability  of  performing  parallel 
tracking  and  identification. 

4.  The  compatibility  with  optical  storage  techniques. 

We  would,  however,  use  the  white-light  or  incoherent  source 
rather  than  the  coherent  optical  counterpart,  because  of  the 
following  advantages: 

1.  It  eliminates  the  coherent  artifact  noise  that 
generally  plagues  the  coherent  optical  system. 

2.  The  white-light  tracking  system  is  generally 
economical  to  operate  and  easy  to  maintain. 

3.  The  white-light  technique  is  very  versatile  and 
simpler  to  handle. 

4.  The  white-light  system  is  generally  not  sensitive 
to  vibrational  problem. 

5.  The  white-light  tracking  system  is  relatively 
smaller  and  simpler  to  install. 

We  emphasize  again,  the  basic  objective  of  this  proposed 
research  is  to  improve  the  tracking  and  identification  capa¬ 
bilities  so  that  reliable  tracking  signals  can  be  obtained. 

1 . 2  Image  Subtraction  with  Extended  Incoherent  Source 

We  shall  in  this  section  describe  a  source  encoding  tech¬ 
nique  for  incoherent  image  subtraction.  We  note  that  this 
source  encoding  techiique  is  capable  of  reducing  the  spatial 
coherent  requirement  of  the  light  source  for'  image  subtraction 
and  at  the  same  time  suppresses  the  coherent  artifact  noise  of 
the  processing  system.  In  other  words,  the  incoherent  image 
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subtraction  system  is  capable  of  subtracting  the  images  in 
complex  amplitude  like  a  coherent  processor,  and  at  the  same 
time  it  suppresses  the  coherent  artifact  noise  like  an 
incoherent  processor. 

Since  the  image  subtraction  problem  is  essentially  a 
one-dimensional  processing  problem,  a  very  narrow  slit  inco¬ 
herent  light  source  can  be  utilized.  Since  the  spatial  coherent 
requirement  for  image  subtraction  is  only  required  at  every 
corresponding  image  points  between  the  two  images,  therefore 
only  a  point-pair  spatial  coherence  between  the  two  images  is 
required.  In  order  to  obtain  a  point-pair  spatial  coherence 
over  the  input  plane  P 2  of  an  incoherent  processor  of  Fig.  1-1, 
we  would  encode  the  extended  source  with  a  set  of  narrow  slit 

apertures  over  the  source  plane  P^.  Now  with  reference  to 

2 

Fig.  1-1  and  the  use  of  the  Van  C i t t er t -Z er nike ' s  theory  ,  the 

degree  of  spatial  coherent  over  the  input  plane  P,  of  the 

3-4 

incoherent  optical  processor  can  be  shown 


(  x  L  -  x  ,  )  = 


x  -x 

S  in  (  ZNt  -  > 


x  -x 
N  S  in  (  2-  ) 


Sin  2  t 


S  (xi'x2^ 

h  d 
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2  h 


2- 


s  (x1  -x? ) 
h~d 


(1-1) 


o  o 

where  d  is  the  spacing  of  the  slits,  s  is  the  slit  width,  and 
N  is  the  total  number  of  slits  over  the  extended  source.  From 
the  above  equation  we  see  that  the  degree  of  spatial  coherence 
for  each  pair  of  points  between  thj  two  images  over  the  input 
plane  P,  separated  by  a  distance  2hQ  can  be  shown  eauals  to 

j  (2h  )  =  Sinc(-~ )  . 
o  d 


( l-:) 
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g 

Thus  for  < <  1 ,  a  high  degree  of  spatial  coherence  between 
the  cori jsponding  image  points  can  be  obtained,  and  a  brighter 
subtracted  image  can  be  observed  at  the  output  plane  of 
the  processor.  We  see  that  the  incoherent  processor  is  indeed 
capable  of  processing  the  image  in  complex  amplitude  like  a 
coherent  processor.  In  comparison  with  a  single  slit  case^, 
the  intensity  of  the  output  image  in  increased  N  fold,  so 
that  the  output  subtracted  image  can  be  made  observable.  We 
would  note  that  the  basic  objective  of  source  encoding  is 
to  provide  an  appropriate  coherent  requirement  for  certain 
optical  information  processing  operation  such  that  more 
efficient  utilization  of  source  power  can  be  achieved. 

We  would  also  note  that  the  source  encoding  technique  can 
be  extended  for  extended  white-light  source,  a  research  program 
is  currently  still  pursuing.  In  addition,  pseudocolor  encoding 
technique  can  also  be  implemented  for  image  subtraction  for 
target  tracking  position.  In  other  words,  it  is  possible  to 
deplay  the  subtracted  image  (i.e.,  target)  in  one  color  super¬ 
imposed  with  the  background  in  another  color. 

Although  the  proposed  incoherent  image  subtraction  tech¬ 
nique  is  aiming  at  the  application  to  missile  tracking  and 
identification,  there  are  several  other  areas  of  intense  intere 
to  the  U . S .  Government  relying  on  signal  tracking  technology 
and  hence  it  stands  to  benefit  from  our  research  work.  Those 
areas  are  radar,  sonar,  etc. 

1 . 3  Report  Outline 

We  deal  primarily  on  an  image  subtraction  with  incoherent 
source  as  applied  to  tracking  and  identification.  Since  the 
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output  plsne  P  ^  will  bo  searched  for,  that  is 


J(x*,y';x',y*;  A )  =//// J '  (x1 ,  ;  *2 ,  y  ;  A )  exp  -  j— j+y  ’  3  1 


_x2:t2'y232)  }da1dS1da2dS2 


(2-191 


Basic  interest  is  usually  centered  on  the  output  intensity 
distribution  I(x',y';  \);  by  letting  x'i=x'a=x'  and 
y'^=y'2=y'  Eq.(2-19)  can  be  written  as 


;  <x  =//v  (x0 -yo;x)  dxodyo-^s  ( V’V-'-  f,J0  ’  yo+i  i }  K  ( v £  i } 

x  expi-j^^x'H^y’ )  }daidli;/S*  (-x+a  .yo+i,) 


x  H  (o2,32)exp( jfj(a2x,+32y,)}da2d22 


(2-20) 


3y  interchanging  the  variables  ( a  ^  ^  )  and  ( a t  t  "  ?  )  Co 

(  i , ~  )  ,  Eq.(2-20)  reduces  to 


!(:■:'  ;  •.)=//y(x  ,y  ;  \)  /;S  (x  r,  ,y  +£ )  H(a ,  £ ) 

'  '  n  n  ■  '  ‘  Q  JO 


exp-.  -  ifc(ix 1  +3y ' )  •  dadc  ~dx  dv 
\  r  o  '  o 


(2-21) 


We  define  S(')  as  the  relative  spectral  intensity  of  the 
light  source  and  C()  as  the  relative  'spectral  response 
sensitivity  of  the  detector  or  recording  material.  Since 
the  different  wavelengths  of  the  light  source  are 
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J(a  ,3  ;ut  ,81;X)=/j’y(x  ,y  ;.X)lS(x  -Hi  y  -hi  )+ -^S(x  -Hi  -Xfv  ,y  +3  ) 
±122  "  oo  o  1  o  1  2  o  1  o  o  1 


+  is(x  -Hi  HXfv  y  +3  ) ] [ S*  (x  Ha  y  HS  )H-k*(x  Ha 
2ol  ool  o2o22  o2 


1  * 


HXfv  ,v  +3  )+  — S  (x  Ha  -Xfv  ,y  H3  ) ] dx  dy 
o  o  2  2  o  2  o  o  2  oo 


(2-15) 


where  the  integration  is  over  the  source  plane  P  ,  and 
S(a,3)  is  the  Fourier  spectrum  of  the  input  s(x,y)  i.e., 


S(ct,3)=F\  s(x,y); 


(2-16) 


If  the  complex  transmittance  of  the  spatial  filter  at  the 
Fourier  plane  is  equal  to  H  0*  ,2),  the  mutual  intensity 
function  immediately  behind  the  spatial  filter  will  be 


J’  (a 


,S2;X)=J(a1 ,31;a2,32;\)H(a.,31)H  (o^.S,) 


(2-17) 


For  convenience,  assume  that  the  spatial  filter  is  only 
effective  over  one  of  the  smeared  Fourier  spectra;  then 
Eq.(2-15)  becomes 


J’  (a 


l’-l 


,;X)=jjy(x  ,y  ;.X)S(x  Ha  -Xf\ 
:  ' J  o  o  o  1 


:  )S  (x  Ha  -Xfv  , 
1  o  2  o 


V 

'  o 


H£1)H(a1,£1)K 


(a 


2  ’  “2 


,)dx  dv 
o  -  o 


(2-18) 


where  the  appropriate  constant  has  been  ignored. 

In  the  final  step,  the  mutual  intensity  function  a 


For  simplicity,  assume  that  the  achromatic  lenses  are 


infinitely  extended.  By  Eq.(2-ll),  the  mutual  intensity 
function  at  the  input  plane  due  to  the  sojrce  intensity 
distribution  ^ (x0,yo)  for  every  wavelength  can  be  written 

J(x1,y1;x2,y2;X)=//Y(xo,yo;X)exp{-j^j1[  (x^-x^x^yj-y^y^  }dxQdyo  . 

(2-12) 

where  (x^,y^)  and  (  *2  '  ^  2  ^  are  cwo  arbitrary  points  over 
the  input  plane  (x,y),  f  is  the  focal  length  of  the 
achromatic  transform  lens,  and  the  integration  is  over  the 
source  plane  PQ.  The  mutual  intensity  function 

immediately  behind  the  diffraction  grating  G  can  be 
written 

J’ (x1yi;x2,y,;X)=J(x1,y1;x2>y,;X)s(x1,yi)s*(x,,yi) 

x(  1+cos  (2-tv  x  )  ]  [l+cos(2rrv  x  )  ]  ,  (2-13) 

o  i  o  2 


where  the  superscript  denotes  the  complex  conjugate,  and 
J  is  the  spatial  frequency  of  the  diffraction  grating. 
Similarly,  the  mutual  intensity  function  at  the  Fourier 
plane  ?2  can  be  written 

J(a1,31;a1,32;X)=j )  <  j  J'  ^ , y  ;x2  ,y2 ;X ) exp{  - x.+£  y , 


-a,x2-32y9]}dx1dyidx2dv0  , 


(2-li) 


where  the  integration  is  over  the  input  plane  P^.  Sv 
substituting  Eqs.(2-12)  and  (2-13)  into  E  q  .  (  2  -  1  4  )  ,  we  have 
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Consider  the  partially  coherent  optical  processing 
system,  which  was  proposed  by  Yu  diagrammed  in 

Fig. 1-3.  Due  co  the  diffraction  grating  the  spectrum  of 
the  input  signal  will  disperse  into  rainbow  colors  in  the 
Fourier  plane  P2  achieving  higher  temporal  coherence  along 
the  a- direction.  Therefore,  the  input  signal  s ( x , y )  can 
be  processed  in  complex  amplitude.  If  the  processing 

operation  is  one-dimensional,  a  fan-shaped  spatial  filter, 
which  is  shown  in  Fig. 2-2,  can  be  used  to  accommodate  the 
wavelength  variation.  For  a  two-dimensional  processing 

operation  a  set  of  narrow  strip  filters  H  n  ^  a  n  •  ~n  ^  1 

(n-1,2, . n)  ,  can  be  used  to  process  different  color 

signals  band  by  band  as  shown  in  Fig. 2-3  Each  of  the 
spatial  filters  is  synthesized  for  a  narrow  spectral  band, 
for  which  the  center  wavelength  is  \ Qn  and  the  bandwidth 
is  equal  to  1'  n.  The  filtered  signal  will  be  Incoherently 
superimposed  at  output  plane  P^. 

The  calculation  of  the  intensity  distribution  at  the 
output  plane  may  be  divided  into  three  steps  as  follows 

(1)  Derive  tne  output  mutual  intensity  function  for  a 

single  wavelength  '  . 

(2)  Calculate  the  output  intensity  distribution 
In(x',y'),  which  is  related  to  the  narrow  strip  spjtiji 
filter. 

(3)  Add  I(x',y')  (n  through  l  to  N)  to  obtain  the  total 
output  intensity  distribution  i ( x '  , y ' )  . 
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where  Y  (x0,yQ)  is  the  intensity  distribution  of  the 
incoherent  light  source,  and  ( xQ , y 0 )  is  the  coordinate  of 
the  light  source  plane.  In  this  case  the  mutual 
intensity  function  at  the  input  plane  is  simply  equal  to 


j  (x  -x2 ; y1-y2)-//y (x  , yQ ) exP{- ( Vx2)xo+(yry2)yo ] } dx0dyo  • 

oo 

(2-11) 

It  is  evident  from  Eq.(2-ll)  that  the  mutual  intensity 
function  at  the  input  plane,  illuminated  by  an  extended 
incoherent  source,  is  a  space  invariant  function. 


2  .  3  General  Formula  for  a  Partially  Coherent  Optical 
P  r  o  c  e s  sor 


An  optical  information  processing  operation  under  the 
partially  coherent  regime  will  now  be  described.  The 

optical  processing  operation  depends  on  the  degree  of 
spatial  and  temporal  coherence  for  the  source,  as  do  the 
spatial  frequency  response  (  i . e . , t r an s f e r  function)  and 
the  noise  performence  of  the  optical  processor.  We  will 
investigate  the  coherence  requirements,  of  a  general  and 
specific  optical  processor.  We  will  look  into  the  e*  :  ects 
of  the  size  and  the  sDectral  bandwidth  of  the  light  source 
e  m  p  loved  on  the  intensitv  distribution  at  the  output  plane 
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Thus  the  mutual  Intensity  function  of  the  light  at  the 
print  pair  <  a  .  3x  5  and  2*^2^  is  rcloted  to  the  Fourier 
spectrum  of  the  mutual  intensity  at  the  input  plane.  We 
note  that  the  Fourier  transform  relation  between  the  input 
and  the  focal  plane  for  the  mutual  intensity  function  is 
not  an  exact  one,  due  to  the  presence  of  the  quadratic 
phase  factor  that  precedes  the  four-dimensional  integral. 

( b)  If  we  let  dQ  be  equal  to  f,  the  quadratic  phase 
factor  vanishes  and  Eq.(2-8)  could  be  again  simplified  as 


X 

J.(a  ,3  ;o  3  )<//// 

±  GO 


Jl(xl,yi'x2,y^exp  ^2?^  (alxl+3lyl) 


-(a,x  +3^)]  )dx  dy  dx^dy^ 


(2-9) 


Evidently  when  the  input  is  placed  in  the  front  focal 
plane  of  an  ideal  thin  lens,  an  exact  four-dimensional 
Fourier  transform  relation  can  be  obtained  for  the 
mutual  intensity  function. 

(c)  In  partially  coherent  optical  processing  an  extended 
incoherent  light  source  is  usually  placed  at  the  front 
focal  plane  of  the  collimation  lens  to  illuminate  the 
input  transparency,  which  is  inserted  in  the  back  focal 
plane  of  the  collimator.  Then 


3 Otr  . > *2 » y 2^  =  i 


,  y(x  ,y  ) 
j  o  o 


for  x  =x  =x 
1  2  o 


otherwise 


(2-10) 
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wher  c  Is  a  constant, 

t<|.(2-7)  Illustrates  that  the  mutual  Intensity  of  the 
light  at  the  output  plane  could  be  calculated  by  a  four- 
dimensional  integration. 

To  this  point  we  have  entirely  neglected  the  finite 
extent  of  the  lens  aperture.  We  note  that  such  an 
approximation  is  an  accurate  one  if  the  distance  dQ  is 
sufficiently  small  to  place  the  input  deep  within  the 
region  of  Fresnel  diffraction  for  the  lens  aperture.  This 
condition  is  satisfied  in  the  vast  majority  of  problems  of 
interest.  The  limitation  of  the  effective  input  by  the 
finite  lens  aperture  is  known  as  a  'vignetting'  effect. 
Note  that  vignetting  in  the  input  space  is  minimized  when 
the  input  is  placed  close  to  the  lens  and  when  the  lens 
aperture  is  much  larger  than  the  input.  In  practice  it  is 
often  preferred  to  place  the  object  directly  against  the 
lens  in  order  to  minimize  vignetting. 

There  are  three  particular  cases  of  the  general 
formula  (2-7)  worth  mentioning: 

(a)  For  the  case  that  d^»f,  the  output  is  placed  at  the 
back  focal  plane.  Eq.(2-7)  becomes 


‘J4^al’^l’a2’^2^=C 


■» 

Jj/JVVVV 


exp{  j^y-(l-  -jr)  [  (ai+-p"(a2+32^  ^ '' 
y^expf-jljKo^+a^J-Co^+a^)]:.  lx1dy1dx2dv 


X 
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°2  .  X2 


2222  ai  xi  Bl  yi 

x(cifnr52  it  +  t  d7 +  <r  )+2?2(  17 +  d- ) 

1  o  1  o  1  o 

G2  ,  y2 


+2n2(  —  +  —  )]}d51dn1dF,2dn2]dx1dy1dx2dy2 
al  o 


(2-5) 


note  that 


Ct  X 

j  exp{  jk[  (  7-  +  j-  -  7  )«?-25,(  7^  +  )  ]  d5,=exp{-jx(  -r±  +  ^  ) 

-«  d1  dQ  f  1  1  d1  dQ  1  d1  dQ 


“l  .  X1 


'<  7-  *  7-  -  7  >  +  r  -  f  >tEr<  17 +  r  >/(  37  -  r 

1  0  1  o  1  o 


d  '  d  f 
1  o 


JkVr- f 

1  o 


f  )]2}dCi“  i[  „;,i  'A  ""1;  ]2exp{-jk(  ^ +  ^ )2/(  ^ +  d; '  i } 


(2-6) 

There  are  similar  formulas  for  variables  n^,  and  1  j ■ 
Therefore  the  transformation  formula  for  the  mutual 
intensity  function  due  to  an  ideal  thin  lens  is 


Vai’WS2)=C  exp{j  +  1  )][(aj+3j)-(a‘+3‘)]; 

1  1  dj_  dQ  f 

x////  J1(x1,y1;x2,y2)exp{j  -^“[l - j — j — —  ]  [  (x“+y") - (x'+yp  ]  -l 

OD  "  o  i  ' - ‘ —  — ' ' 

(-lWl)~(;2W2) 


9  9  9  9 


o  d 

0  d,  d  f 
1  o 


x  exp’-jk - — f— - — ; — — — — —  }dx.  dv.  dx.dy? 

d  d,(  7-  +  7-  -  7  ) 

old.  d  f 
1  o 


(2-7) 
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dimensional  Fourier  transformation  of  a  light  signal  in  complex 

14  15 

amplitude  ’  .  For  partially  coherent  light  we  can  show  that 

the  four-dimensional  Fourier  transformation  of  the  mutual 

intensity  function  could  be  performed  by  an  ideal  thin  lens. 

It  can  be  shown  that  a  plane  wave  passing  through  a  thin 

lens  experiences  a  phase  delav  factor  of  t  (5,n)  such  as^, 

o 

to(5,n)  =  exp[jknAo]expMkU2+n2/2f]  ,  (2-1) 

where  n  is  the  refractive  index  of  the  lens  material,  f  is  the 
focal  length,  Aq  is  the  center  thickness  of  the  lens,  k  =  2  tt  /  X  , 
where  A  is  the  wavelength  of  light,  and  (5,n)  is  the  coordinate 
system  of  the  lens  plane.  The  first  term  of  Eq .  (2-1)  is  simply 

a  constant  phase  delay,  while  the  second  term  may  be  recognized 
as  a  quadratic  phase  factor  associated  with  a  spherical  wave. 

Consider  the  more  general  geometry  of  Fig.  2-1.  where  an 
object,  located  a  distance  dQ  in  front  of  the  lens,  is  illumin¬ 
ated  by  spatially  partial  coherent  light.  If  the  mutual  intensity 
function  at  the  object  plane  is  J ^ (x^ , y^ ; x^ , Y ^ ) ,  our  goal  is  to 
try  to  find  the  mutual  Intensity  function  at  the  output  plane, 
which  is  behind  the  lens  at  distance  d^.  The  light  at  the  lens 
plane  (  ; ,  n )  can  be  expressed  by  using  the  superposition  integral 
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of  correlation  in  an  investigation  relating  to  image  formation 
for  •  microscope^. 

An  important  stage  in  the  development  of  this  subject 

9  10 

was  made  by  Van  Cittert  and  Zernike  .  They  determined  the 

"degree  of  coherence"  for  the  light  disturbances  at  any  two 

points  on  a  screen  illuminated  by  an  extended  primary  source. 

The  method  was  simplified  further  and  applied  to  the  study 

of  image  formation  and  resolving  power  by  Hopkins^.  Wolf 

introduced  a  correlation  function,  which  was  found  to  obey 

rigorously  the  wave  equation,  to  depict  the  correlation  proper- 

12 

ties  of  the  light  disturbances  at  any  two  points  .  Such 
'second  order'  correlation  functions  are  adequate  for  the 
analysis  of  the  usual  optical  experiments  involving  interference 
and  diffraction  of  light  from  steady  sources. 

Before  investigating  the  behavior  of  a  partially  coherent 
optical  processor,  it  is  necessary  to  establish  the  transforma¬ 
tion  relations  of  the  mutual  intensity  function  (or  correlation 
function),  which  is  equal  to  the  degree  of  coherence  in  the 

image  of  an  extended  incoherent  q ua s i -mo no c hr oma t ic  illumination. 

1  2 

In  this  section,  we  use  Wolf's  theory  of  partially  coherent 

light  to  develop  a  transformation  formula  for  the  mutual  intensity 

function  passing  through  an  ideal  thin  lens.  We  then  apply  it  to 

derive  the  general  formula  for  the  partially  coherent  optical 
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processor  which  was  introduced  by  Yu 

t 

2 . 2  Propagation  of  Mutual  Intensity  Function 

One  of  the  most  remarkable  and  useful  properties  of  a 
converging  lens  is  its  inherent  ability  to  perform  the  two- 


image  subtraction  operation  relies  on  complex  amplitude  opera¬ 
tion,  the  partial  coherent  requirement  of  the  extended  source 
linage  subtractor  should  be  evaluated.  The  feasibility  and  the 
constraints  of  the  svstem  performance  should  be  discussed.  We 
shall  report  in  the  following  the  analysis  of  the  coherence 
requirement  for  the  proposed  incoherent  image  subtraction  system. 
The  application  of  this  proposed  image  subtraction  system  to 
IC  chip  inspection  are  given.  Experimental  demonstrations  for 
the  IC  mask  inspection  is  included.  We  shall  also  report  a 
technique  of  generating  a  2-D  color  hologram  with  the  image 
subtraction  technique.  An  experimental  result  of  color  holo¬ 
gram  image  obtained  with  this  proposed  system  is  also  provided. 

1 1 .  Image  Formation  by  Partially  Coherent  Illumination 
2 . 1  Introduction 

Researchers  in  optical  information  processing  have  confined 
their  attention  to  the  analysis  and  synthesis  of  optical  systems 
in  which  the  illumination  is  either  extremely  coherent  or 
incoherent.  Surely  we  must  expect  a  continuous  transition 
between  these  two  extreme  limits.  Such  a  transitional  region 
exists  and  is  known  as  the  field  of  partial  coherence. 

The  first  investigation  which  has  a  close  bearing  on  the 
subject  of  partial  coherence  was  done  by  Verdet.  who  studied 
the  size  of  the  region  of  coherence  for  light  from  an  extended 
source^.  Later,  Michelson  established  the  connection  between 
the  visibility  of  interference  fringes  and  the  intensity  dis¬ 
tribution  on  the  surface  of  an  extended  primary  source^.  A 
further  contribution  was  made  by  Berek,  who  used  the  concept 
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mutually  incoherent,  the  intensity  distribution  related  to 
th  narrow  strip  filter  can  be  written  as 


X  +AX  / 2 

on  n  30  => 

In(x'>y')=  /  J/y (xQ , yo ; A ) S  (X )C (X ) ! //s (xQ+a-X  f v  , y  +S ) 

^on  °° 


({(a.SjexpC-j-rrCx'a+y' 3)  }dad6  [  ^dx  dy  dX 
a  1  00 


(2-22) 


where  X  is  the  center  wavelength  and  AX  is  the 

on  n 

bandwidth  corresponding  to  the  nth  narrow  strip  filter. 
Finally  the  overall  intensity  distribution,  due  entirely 
to  the  spectral  band  of  the  light  source,  at  the  output 
plane  would  be. 


I(x',y')=V  I(x\y’) 


n-1  “ 

N  Xon+iAn/2  - 

=  I  ^-AX  / 2  ^y(V yo;X)S(X)C(X)  //s<K  4*--. 
n=l  on  n  00  * 


fv  ,  y  -3) 
o  o 


H(a, £)expt -jy-~(x' a+v ' 3) ;dad£ , “ax  dv  i\ 

A  I  "  O  O 


C-23) 


which  is  essentially  the  incoherent  summation  of  the  filtered 
s igna Is . 

III.  Coherence  Requirement 
J  .  1  Introduction 

Since  its  invention,  the  laser  has  been  a  fashionable 


tool  for  many  applications,  particularly  for  coherent  optical 


21 


information  processing.  This  trend  is  mainly  due  to  the 
c  'tnplex  amplitude  processing  capabilities  of  the  coherent 
light.  Unfortunately,  coherent  processing  systems  are 
also  contaminated  with  coherent  artifact  noise,  which  is 
regarded  as  the  number  one  enemy  of  coherent  processings. 

In  addition,  coherent  sources  are  expensive  and  the  pro¬ 
cessing  environment  is  usually  critical. 

Many  optical  information  processing  operations  can 
be  implemented  by  systems  that  use  extended  polychromatic 
light  sources;  however,  the  schemes  are  not  without  pit- 
falls.  The  partially  coherent  processing  system  is 
capable  of  reducing  the  inevitable  artifact  noise,  but 
generally  introduces  a  dc-bias  build-up  problem,  resulting 
in  poor  processing  performance. 

We  have  shown  that  image  subtraction  can  be  achieved 

3  4 

by  the  use  of  an  extended  incoherent  source  ’  .  The  question 

to  be  addressed  is,  to  what  degree  can  we  relax  the 
coherence  requirements  without  sacrificing  the  processing 
capabilities  of  the  processing  system.  The  nature  of  each 
optical  processing  operation  governs  the  temporal  and 
spatial  coherence  requirements  necessary  to  obtain  satis¬ 
factory  results.  We  will  reduce  the  general  formulas  we 
have  obtained,  in  the  previous  section,  to  suit  the  image 
subtraction  application  under  consideration.  In  the 
following  we  will  discuss  the  temporal  and  spatial  coherence 
requirements  for  the  image  subtraction  operation. 
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1 . 2  Coherence  Requirement  for  Image  Subtraction 

In  a  recent  paper  published  by  Wu  and  Yu  [3],  a 
technique  of  image  subtraction  with  an  encoded  incoherent 
source  is  described.  Fig.  3-1  illustrates  the . processing 
system  for  image  subtraction.  An  extended  incoherent 
source  with  an  encoding  mask,  is  placed  at  plane  PQ;  the 
two  object  transparencies  to  be  subtracted  are  located  at 
plane  P^.  A  sinusoidal  grating  (either  a  phase  or 

amplitude  grating)  in  Fourier  plane  ?2  acts  as  a  spatial 
filter  to  perform  the  subtraction  operation.  With  the 
reduced  spatial  coherence,  image  subtraction  can  also  be 
obtained  with  a  source  encoding  mask.  In  this  section  the 
coherence  requirments,  which  depend  on  the  source  size  and 
the  spectral  bandwidth  of  the  lighc  source,  will  be 
evaluated  for  the  image  subtraction. 


j.2.1  Basic  Formulas 

The  approch  will  consist  of  calculating  the 
propagation  of  the  mutual  intensity  function  through  the 
optical  processor  chat  utilizes  an  encoded  source  of 
extended  size.  With  reference  to  the  partially  coherent 

I 

optical  processing  system  of  Fig.  3-1,  the  encoded  source 
intensity  distribution  can  be  written 
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Y ( x  y  )=rect(  — ,  )’ 
o  o  a 


N 

'  7  5(x  -nD) 
n— N  0 


(3-1) 


where  *  denotes  the  convolution  operation,  d’  is  the 
spacing  of  the  encoding  slits,  d'  is  the  slit  width  and 
2S+1  is  the  total  number  of  encoding  slits,  as  depicted  in 
Fig.  3-2.  Since  the  image  subtraction  is  a  o ne- dimen s ion al 
processing  operation,  a  on e- dimen sio nal  notation  will  be 
adopted  for  the  discussion.  Thus,  from  Eq.(2-22),  Che 
encoded  mutual  intensity  function  at  the  input  plane  can 
be  written  as 


®  x  N1 

J(x  — x,;\)=/  rect(  )*  l  5  (x^nDWp  [-jfj(x1-x2)*Jdx0 
1  -®  n=-N 


=  /rect  (  )exp[-jIj(x1-pc2)xoidxo 


x  /  l  5(xo-nD)expi-j~-(x1-x2)xo]dxo 
- »n=-N 


N  dUj-xJ 

*  i  sine  — 5jy^-exp[-jYj(x1-x2)nD] 
n=-N 


(3-2) 


The  two  object  transparencies,  in  Fig.  3-1,  A(x)  and  B(x) 
are  placed  at  the  input  plane  They  can  be  expressed  as 
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s (x)=A(x-h)+B(x+h)  , 


(3-3) 


where  2h  is  Che  separation  between  the  transparencies. 
From  Eq.  (3-2)  the  mutual  intensity  function  at  the 
spatial  frequency  plane  Pj,  is 


M  co  ' 

,  d (x  *x  ) 

J(i. ,a„;X)=£  //  sine [- 

1  -  n=-M 


Xf 


-]exp[-j|j(x1-x2)nD] 


X 


[A(x^-h)+B  (x^+h)  ]  [A  (x2“M+3  (x^+h)  ] 


x 


'  T 

expC-j^j  (a1x1-a2x2) ]dxxdx2 


(3-1) 


In  the  image  subtraction  operation,  a  sinusoidal  gracing 
with  spatial  frequency.)  is  placed  at  the  Fourier  plane. 
The  mutual  intensity  function  immediately  behind  Che 
grating  is 


(a, ,a,;X)=J(a, ,a.;i ) (l-csin2nv  a  ) ( l+esin2 
12  12  o  I 


’v  a,) 
o  2 


(3-5) 


The  intensity  distribution  at  the  output  plane  P-j,  for  a 
given  *  Is  therefore 


Kx’  ;  >.  )=  ;  j'  (a  ,a,  ;\  )exp[-j^j(a  -a2)x'  Idc^dc^ 


(  3-6  ) 


By  using  Eqs.  (3-H  and  (3-51,  Eq.  (3-6)  becomes 
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l(x’;0  =  l  / /sine  [ - j^-p-]exp[j— (Xj-x^nD] 

n=-N 

x  [A(x^-h)+B(x^+li)  ]  [A  (x2~h)+B  (x^+h)] 

x  //  (l+csin2irvoa  )  (l+csin2irvoa2)exp  [-j-^a^x^-c^x,,)  ] 
x  exp  [-jr^a^-a^x'  ]da2da2}dXidx2 


Ic  is  easy  to  show  that  the  second  integral  in  Eq.  (  3-7) 
is  equal  to 


2- 

ff (l+csin2Tu  a  ) (l+csin2iv  a  )exp[-jf-r(a, x, -a_x  ) ] 
oi  o  2  3 r  l 


X  exp[-j-|y(a  -a2)x' ]da1da2 


-[5(xi-x,)+  2^3(x1-x'-\fvo)-  I^(x2-x’+\fvo)] 


[  £  ( x,-x  *  )+2^  (x,-x ’  -  \  f vq  )  -  2f:  (x2-x  ’  +\  f  vo )  ] 


(3-8) 


II y  substituting  Eq.  (3-3)  into  Eq.  (3-/),  and  after  a 
tedious  calculation  I(x";  \) 


I(x'  ;\  )=>!(  A(x'-h)  |+  3(x'+h)|  “) 


can  be  written 
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-*•  “[2sinc  (2ndv  )Re{A(x' -h+Xfv  )B(x'+h-Xfv  ) 

4  0  0  O 

-;A(x* -h+Xfv  ) |2-!B(x'+h-Xfv  )|2 
1  o  o 

-!A(x'-h-Xfvo)|2-iB(x'+h+X£vo)|2  ,  (3 

where  Re{f(x)}  represents  the  real  part  of  function  f(x). 
Eq.  (3-9)  illustrates  that  there  will  be  six  difraction 
terms  at  the  output  plane  if  we  choose  u  * h/ X f . 
Considering  only  the  diffraction  terms  which  are  around 
the  optical  axis,  we  have 

1^  (x'  ;X)  =  j  A(x'  -h+Xfv^)  I  2-2sinc (2ndvQ) 


x  Re •.  A(x ' -n+X f v  )B(x'+h-Xfv  ) 
o  o 


+■ B(x'+h-Xfv  )!2 
1  o 


(  3-10) 


From  Eq  .  (3-10)  it  is  noted  that  if  the  slit  size  d  equals 
zero,  the  analysis  reduces  to  the  case  of  strict  spatial 
coherence  with  the  intensity  distribution'  given  by 


I(0)(x';\>:dl0  =  [ A ( x ' )  - 3 ( :< ' )  j ' 


{  3-11) 


Eqs.(3-10)  anti  (3-11)  show  that  a  perfect  subtracted  image 


requires  a  strictly  spatial  coherent  system,  and  the 
quality  of  the  subtracted  image  decreases  as  the  slit  size 
increases. 

To  analyze  the  case  of  partial  coherence  assume  that  the 
light  source  has  a  uniform  spatial  bandwidth  and  that  the 
spectral  response  of  the  detector  is  also  uniform 
[i.e.,S(X)»k;C(X)«k].  In  this  case  the  image  intensity 
distribution  at  the  output  plane  may  be  given  by 


I(0)(x> 


Xo+AX/2 

J  K* 

X  -AX/2 
o 


;  x  )d  x 


X o+AX / 2 

=  /  !A(x'-h+Xfv  )  1 2-2sinc(2rrdv  ) 

‘  o  o 


x  Re { A ( x ' -h+X f v  )B(x'+h-Xfv  )} 
o  o 

,  a 

+'3(x'+h-Xfv  )  "  dX 
o 


H-12  ) 


For  polychromatic  ilium  in ation  the  spatial  frequency  of 
the  diffraction  grating,  v  ,  is  chosen  to  be  uQ»h/Xf, 

where  \  is  the  center  wavelength,  and  AX  is  the  spectral 

a 

bandwidth  of  the  light  source.  Changing. the  variable  X  to 
\  '  m  X  ,  Eq  .  (  3  -  1 2  )  becomes 
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AX/2 

I  ^  (x'  )=  /  {[  A(x'-X'fv  )  J  “-2sinc  (2ti<Jv  ) 

o  o 

-AX/2 


x  Re[A(x'+X'fv  )B(x'-X'f\)  )] 
o  o 

+  lB(x’-X’fv  )  1 2>dX '  .  (3-13) 


Eq.  (3-13)  may  be  simplified  by  using  a  Taylor  series 
expansion  for  the  input  object  functions,  i  .  e  .  > 


A(x'+X'£v  )=A(x')+  ) 

0  m=l 


l  JL  (^)  /  ,  \  |  -  \  HI 

—  A  (x  )  (X  fv  ) 
m!  o 


(3-14a) 


and 


B(x' -X ' fv  )=B(x')+  " 

0  “  i 

m=l 


(-1)* 

m! 


B(m) 


(x'HX'fv  )m 
o 


(  3  —  1 4  b  ) 


where 

.  ( m  .J  .  r  \  ^  ^  i  \  /  .  , 

A  (x  )  =a  x  )  /  ax 


and 


i  ^  (x'  )=dTn3 (x* ) /dx,:a 


Therefore  we  find  that 


I(0;(x')*‘ .A(x') ;2+'sinc(2^dv  )A(x' )B(x' )*' 3(x' ) | 2  h 


(  3  - 1  5  a  ) 


(  3  -  1  5  b  ) 
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m=even 

m=0 


2  (m+1 ) ! 


ffv  )m[A^m^  (x '  )  -sine  (2?idv  )f/m^(x')] 


x  [A(x')+B(x’) ] (AA ) 


m+1 


mm'  2  m  (m+m’+l)m!m' ! 

m+m' =even 
m=0 
m'=0 


(fv  )nWm'{A(m)(x')A(m')(x') 
o 


-(-l)m’siRc(2-Tdv  )A(m)  (x')B(n,)(x') 
o 

+(-Drt’B(ai')(x')}w)rt,+i  . 


(3-16) 


This  equation  shows  that  a  high  contrast  subtracted  image 
can  be  obtained  with  object  transparencies  of  moderately 
low  spatial  frequency  content.  In  addition  this  equation 
may  be  used  to  compute  the  spectral  requirement  of  the 
light  source. 

An  example  analysis  will  now  be  presented  to  develop 
the  modulation  transfer  function  ( >!  T  F )  .  These  will  be 
used  to  determine  the  temporal  and  spatial  coherence 
requirements  of  the  image  subtraction  process.  Assume 

that  the  input  object  transparencies  are 

A  ( x )  =  1  ,  (  3  -  1  7  a  ^ 

3  (,  x )  -  ~  i  i*C  cosi.  2".x)  ,  p  1  ^  "  *  '  "  ' 

4-  O 


and 
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»nere  vis  the  spatial  frequency  of  the  test  object.  Then 
the  intensity  of  the  subtracted  image  is 

ax/2  c 

I^(x')=  J  { l-2sinc  (2rav  )(  -q-  +  — °  cos2nv  (x  1  -X  '  fv  )] 

'  o  2.  2.  o 

-AX/ 2 

C  , 

+  [  —  +  — ' 3  cos2tv  (x  ' -X  '  f'JQ)  ]  "}  dX  ' 

i 

C“ 

»[  7-  -  — -  -sinc(2-rdu  )]AA+C  [  rr  -sinc(2~dv  )] 

4  8  o  o  1  o 


x  sinc(rfo  4X  j  AXcos  (2-i'jx’ ) 
o 

C 

+  sinc(2~fv  AX )AXcos2~ (2v)x' 
o  o 


For  this  equation  it  is  evident  that  in  addition  to  the 
basic  frequency  term  [cos(2  :vx')]  ,  there  is  a  second 
harmonic  tern  [  c  os  2  (  2 ,  )  X '  ]  that  appears  at  the  output 
plane.  The  corresponding  contrast  measures  can  be 

respectively, written  as 


2C  [ l-2sine(2dv  )]sinc(w  AX) 
o  o  o _ 


*4 


C2 

— °  -sinc(2dv  ) 

3  o 


(  3 


and 
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v  (2v)  = 
P 


•q  sin(2fvv  i,\ ) 
3 _  o 


Co 


sinc(2dv  ) 
o 


(  3  -  1 9  b  ) 


From  these  equations 

,  it  can  be  seen 

that 

the 

subtracted 

image 

depends 

o  a 

the  spatial  frequency 

o  f 

the  input 

object, 

,  v  ,  and 

the 

spectral  bandwidth 

of  the  1 

i  g  ht  source 

2 A  .  Thus,  a  high  spatial  frequency  input  object  and  a 
broad  spectral  band  of  light  source  will  result  in  a  lower 
contrast  of  the  subtracted  image.  In  contrast,  if  i>  =C 
and  d“0  (i.e.,  a  coherence  case),  the  subtracted  iaage 
produces  a  contrast  reversed  iaage,  The  corresponding 
intensity  distribution  is 


1  (x')=[A(x’)-3(x’)j- 
o 


]  r  *> 

:[i-  -q-  -  -  cos  (2tvx' )  ]“ 


I  _  Co 
4  3  2 


cos(2rvx')-  ^-0cos2-t  (2v)x' 


0-201 


Thus,  the  corresponding  contrasts  of  the  basic  ana  second 
haraonic  frequencies  are 


{•j  )  =4C  /(2-C~) 
:  o  o 


<  j  -  2  i  a  / 


<  (u)=  C  /  (2-C  )  , 

c  no 


3  -  2  1  b  ) 


respectively. 

The  MTr  is  defined  as  the  ratio  of  the  contrast  of 
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mask  (i.e.,d/D).  To  obtain  a  high  degree  of  spatial 
coherence  for  a  subtraction  operation,  a  lower  ratio  of 
d/D  should  be  maintained.  The  dependence  of  the 

t  t 

modulation  transfer  function  on  d/D  is  plotted  in 
Fig.  3-3.  From  this  figure  it  can  be  seen  that  when 

t  i 

d/D**0.3,  the  MTF  approches  zero,  and  the  subtraction 
operation  ceases  to  function.  If  the  ratio  of  d/D  is 

sufficiently  small,  a  higher  .MTF  can  be  obtained.  A  few 
numerical  examples  are  provided  in  Table  3-J. 

From  this  table  it  is  clear  that  for  d’/D'=0.05.  a 
relative!';  high  M77  =  0.33  can  be  obtained  for  the  subtraction 
operation.  To  obtain  a  higher  contrast  of  the  subtracted 
image  a  smaller  ratio  of  d’/D’  is  required. 
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In  addition,  the  MTF  decreases  rather  rapidly  as  the 

t 

r  tio  of  the  slit  width,  d,  to  the  focal  length  of  the 
transform  lens,  f,  (i.o.,d/f)  increases.  Thus,  if  the 
focal  lengths  of  the  collimating  and  transfer  lenses  are 
set  equal  to  one  unit  of  the  horizontal  axis,  Fig.  3-7 
will  represent  the  slit  width  of  the  encoding  mask. 

Since  the  visibility  of  the  subtracted  image  depends 
on  the  modulation  transfer  function,  with  required  MTF  and 
h,  the  appropriate  source  size  can  be  determined.  The 
effects  between  the  MTF,  the  separation  h,  and  the  source 
size  are  tabulated  in  Table  3 --2,  where  the  focal  length  of 
the  Fourier  transform  lens  was  assumed  to  be  300  mm 

Table  3-2  shows  that  a  very  narrow  source  size,  d,  is 
required  for  a  high  modulated  transfer  function.  However, 
such  a  small  incoherent  source  is  difficult  to  achieve  in 

practice.  This  difficulty  can  be  alleviated  with  a  source 
encoding  technique.  The  discussion  of  this  technique 
follows . 


A  multislit  source  encoding  mask  is  used  for 
illustration.  Note  that  the  spatial  period  of  the 
encoding  mask  should  be  precisely  equal  to  that  of  the 
diffraction  grating,  G  (i.e.,  d  =  1 / \> Q ) •  In  other  words,  the 
modulation  transfer  function  of  Eqs.  (3---5,  and  ;3--6)  has 
sinc(2  did).  That  is,  the  spatial  coherence  requirement 
Is  independent  of  the  slit  size,  but  it  depends  on  the 
ratio  of  slit  width  to  the  spatial  period  of  the  encoding 


Figure  3-6.  Apparent  Modulation  Transfer  Function  vs 
the  Separation  h. 
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3.2.3  Spatial  Coherence  Requirement 

Consider  the  case  of  perfect  temporal  coherence 
(l.e.,  JU-0  )  and  spatially  partial  coherence,  where 
Eqs  .  (  3-6  8a)  and  (3-68b)are  of  the  fora 


MTFr 


(2-C“) [l-2sinc(2dv  )] 
_ o _ o 

10-C~-8sinc (2dv  ) 
o  o 


(3-25) 


2-C 


mtf2= 


10-C^-8sinc (2dv  ) 
o  o 


( 3-26) 


Note  that  the  nodulation  transfer  functions  are 
independent  of  the  spatial  frequency  of  the  input  object. 

t 

The  above  equations,  however,  depend  on  the  slit  width  d. 
This  requires  that  the  grating  be  precisely  designed  to 
match  the  separation  of  the  object  transparencies,  i.e. 

v  =h/\  £  . 

o  o 

The  plots  of  the  MTF  versus  the  separation  h,  for 
Eqs. (3-25;  and  (3-26)  are  shown  in  Fig. 3-6.  It  is  clear 

that  the  MTF  is  a  monotonic  function  of  h.  To  obtain  a 

high  contrast  subtracted  image,  the  separation  h  must  be 
reduced.  However,  decreasing  the  separation  between  the 
object  transparencies  limits  the  site  of  input  objects  to 
be  processed. 


us  Values  of  Separatlo 
Separation  Between  the 
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the  Spectral  Bandwidth  of  the  Light 
AX  for  Different  Minimum  Desirable 


Apparent  Modulation  Transfer  Function  for 
Partially  Coherent  Image  Subtraction: 

(a)  Basic  Frequencv: 
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The  noraalized  MTF  curves  of  the  basic  and  harmonic 

freq  encies  are  shown  in  Fig.  3-3.  it  is  obvious  that  the 

contrast  of  the  subtracted  image  decreases  mono t on ic all y 

as  a  function  of  the  object  spatial  frequency.  However, 

che  MTF  of  the  subtracted  image  decreases  as  the  spectral 

bandwidth  of  the  light  source  Increases.  In  other  words, 

the  quality  of  the  subtracted  image  improves  as  the 

spectral  bandwidth  of  the  light  source  and  the  spatial 

frequency  of  the  object  decrease.  We  define  y  to  be  the 

cut-off  spatial  frequency  when  the  MTF  decreases  to  a 

minimuQ  value  C  ,  as  shown  in  Fig.  3-3.  The  value  C 
u  n 

depends  on  the  maximum  resolution  of  the  output  recording 
material.  Fig.  3-4  shows  tne  functional  relationship  of 
the  cut-off  frequency  v  c  and  the  spectral  width  i.\  for 
various  values  of  Cm .  It  is  possible  to  determine  the 
spectral  bandwidth  requirement  from  this  figure.  The 

relationship  between  the  cut-off  frequency  y  ,  the 


spectral  width  — '  and  the  separation  between  the  two  input 

transparencies  h  is  shown  in  Fig.  3-5.  Note  that  the 

spectral  bandwidth  required,  for  a  given  cut-off 

frequency,  decreases  with  increasing  separation  h. 

Table  3-1  illustrates  the  dependence  of  dX  on  v  and  h. 

c 

The  focal  length  of  the  Fourier  transform  lens  was 
selected  as  f»300  mm  for  the  calculations.  It  is  clear 

from  the  table  that,  as  the  spatial  frequency  and  object 
separation  increase,  the  spectral  bandwidth  of  the  light 
source  must  decrease. 
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the  Input  and  output  sinusoidal  objects.  From  Eq  .  (3-19) 

and  \.3-21)  these  are  given  by 


MTF (v)= 


[l-2sinc (2dvA) ] (2-dC“)sinc(fvv  AX) 
_ u _ o _ o 

10-C“-8sinc (2dv  ) 
o  o 


( 3-22a ) 


MTF  (2'j)  = 


(2-C^)sinc(2fvv  AX) 
_ o _ o 

10-C^-8sinc (2du  ) 
o  o 


(3-22b) 


As  previously  stated,  Eqs.(3-68a)  and  (3  —  6  8  b )  will 
allow  the  evaluation  of  the  temporal  and  spatial  coherence 
requirements  for  image  subtraction. 


3.2.2~he  Temporal  Coherence  requirement 


The  case  of  strictly  spatial  coherence  will  be 

discussed  first.  This  requires  that  the  slit  size  d  of 
the  encoding  mask  approach  zero.  Eqs  .(  3  -  2  2  a  )  and  (  3-22b)  are 
then  reduced  to 


MTF ( . ) =s inc ( fuv  AX)  , 
o 


MTF(2v)=sir.c(2:vv  AX) 
o 


( 3-2  3  i 


(3-243 


I 
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IV .  Resolution  Limit 

It  is  known  that  the  resolution  limit  is  affected  by  the 
coherence  requirement  and  the  space  bandwidth  product  of  the 
processing  system.  We  shall  now  discuss  in  detail  these  two 
factors. 

4 . 1  Effect  Due  to  Coherence  Requirement 

We  shall  first  discuss  the  effect  on  resolution  limit 
due  to  temporal  coherence  requirement.  We  note  that  for 
the  case  of  strict  spatial  coherence,  the  slit  width  d  of 
the  source  encoding  mask  should  approach  zero.  Thus  the 
MTF  of  Eqs.  (3- 22a)  and  (3-22b)  becomes, 

MTF(u)  =  s  inc  (ph  — )  ,  (4-1) 

'  0 

MTF(2ui)  =  s  inc  (  2ph  ^-1)  ,  (4-2) 

where  d  =  0,  the  normalized  MTF  curves  of  the  fundamental  and 
second  harmonic  frequencies  are  shown  in  Figure  3-3.  It  is 
obvious  that  the  contrast  of  the  subtracted  image  decreases 
monotonicallv  as  a  function  of  the  input  object  spatial 
frequency,  for  a  given  separation.  The  MTF  of  the  subtracted 
image  also  decreases  as  the  spectral  bandwidth  of  the  light 
source.  In  other  words,  the  quality  of  the  subtracted  image 
improves  for  narrower  spectral  bandwidth  of  the  light  source 
and  lower  spatial  frequency  of  the  object. 

Let  us  define  p  ^  be  the  resolution  limit  of  the  sub¬ 
tracted  image,  where  the  MTF  decreases  to  a  minimum  value 

C  ,  as  shown  in  Figure  3-3.  Figure  3-4  shows  the  functional 
m  a  3 

relations  h in  of  the  resolution  limit  p  and  the  snectral 
*  c 
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width  AX  for  various  values  of  C  .  We  shall  now  determine  the 

m 

spectral  bandwidth  requirement  AX  from  this  figure.  The 
relationship  between  the  resolution  limit  p^,  the  spectral 
width  AX  and  the  separation  of  the  input  object  transparencies 
is  depicted  in  Figure  3-5.  We  see  that  the  resolution  limit 
decreases  monotonically  as  the  spectral  bandwidth  and  the 
separation  of  input  objects  increase.  In  order  to  have  a 
feeling  of  magnitude.  Table  4-1  illustrates  the  relation  of 
Pc>  h  and  AX,  where  the  center  wavelength  Xq  of  the  light 

O 

source  is  assumed  to  be  5461A  and  the  HTF(^)  =  0.3.  Thus, 
we  see  that,  to  obtain  a  high  resolution  of  the  subtracted 
image,  a  narrow  spectral  band  of  light  source  is  required. 

We  shall  now  consider  the  spatial  coherence  requirement. 

The  MTF  of  Eqs  (3-22a)  and  ( 3-2  2b )  are; 

(2  -  C2) [1  -  2  sine ( 2dp  ) ] 

MTF  1  =  - — - - 2 -  ,  (4-3) 

10  -  C  “  -  8  sinc(2dp  ) 
o  r  o 

and 

? 

2  -  C“ 
o 

MTF  ,,  =  -  ,  (4-4) 

10  -  C2  -  8  sinc(2dp  ) 
o  o 

where  A '  =  0,  p  ^  =  h / ( X  f ) .  From  these  equations,  we  see  that 

the  MTF's  are  independent  of  the  object's  spatial  bandwidth, 
however,  they  are  functions  of  the  slit  width  d.  The  corres¬ 
ponding  MTF  versus  the  separation  h,  are  plotted  in  Figure  3-6. 
Thus  to  obtain  a  high  contrast  subtraction  image,  the  separa¬ 
tion  h  should  be  small.  But  decreasing  the  separation  also 
limits  the  size  of  input  objects  to  be  observed.  The  relation¬ 
ship  between  the  MTF,  the  object  separation,  and  the  slit 
width  is  tabulated  in  Table  4-2,  where  the  focal  length  of 
the  Fourier  transform  lens  is  assumed  to  be  300  mm. 


From  this  table,  we  see  that,  to  maintain  an  adequately 
high  MTF,  a  very  narrow  source  is  required.  However  a  high 
intensive  narrow  source  is  difficult  to  obtain  in  practice. 
Nevertheless,  this  problem  can  be  alleviated  with  a  source 
encoding  technique,  as  illustrated  in  the  following: 

A  multislit  source  encoding  mask,  is  used  for  this  purpose. 
The  spatial  period  of  the  encoding  mask  is  precisely  equal  to 
the  diffraction  grating  G,  (i.e.,  D  =  1/p).  The  spatial 
coherence  requirement  though  independent  on  the  slit  size, 
is  governed  by  the  ratio  of  the  slit  width  to  the  spatial 
period  of  the  encoding  mask,  i.e.,  d/D.  The  ratio  of  d/D 
is  adequately  small  to  achieve  a  high  degree  of  point-pair 
spatial  coherence  at  the  input  plane.  The  dependence  of  the 
MTF  upon  d/D  is  shown  in  Figure  3-8.  It  is  obvious  that  the 
subtraction  operation  ceases  when  the  MTF  approaches  zero, 
i.e.,  d/D  =  0.3.  A  few  numerical  examples  are  presented  in 
Table  1-3. 

V .  Application  to  Micro-Circuit-Chip  Inspection 
3 . 1  Effect  Due  to  Space  Bandwidth  Product 

Since  the  application  of  image  subtraction  tc  micro-circuit- 
chip  under  inspection  is  a  high  spatial  frequency  type  object, 
the  transform  lenses  L,  and  L are  required  to  be  highlv  dii- 

i  <: 

traction  limited.  Here  are  a  few  commercial  microscopic  objec¬ 
tives  suitable  for  this  purpose,  as  tabulated  in  Table  5-1. 

The  size  requirement  of  the  micro-circuit-chip  under  inspection 

(i.e.,  S  =  h  =  j(4»  2  )  and  the  resolution  limit  of  the  pro¬ 

cessing  system  can  be  calculated  if  a  specific  microscopic 
objective  lens  is  selected.  A  few  numerical  examples  are  given 


in  Table  5-2. 


For  instance,  if  the  spectral  bandwidth  of  the 


incoherent  source  is  10A,  the  slit  size  of  the  source  encoding 
aas  .  is  3  ym  and  the  No.  1  micro-scopic  objective  lens  is 
used,  then  the  optical  processor  has  the  capability  of 
observing  a  micro-circuit  chip  width  as  large  as  9  mm  with 
a  spatial  resolution  as  high  as  220  lines/mm.  The  capability 
is,  of  course,  dependent  upon  the  minimum  detectable  value 
o  f  MTF ,  which  is  less  than  7%  in  the  above  example. 

5 . 2  Experimental  Demonstration 

In  experimental  demonstration,  we  insert  a  standard 
simulated  micro-circuit  chip  transparency  in  one  of  the 
input  apertures  of  the  processor,  as  shown  in  Figure  5-la, 
and  a  defective  or  faulty  one  in  the  other  aperture,  as 
shown  in  Figure  5-lb.  3y  the  comparison  of  these  two  Figures, 
we  see  that  there  are  several  links  missing  in  Figure  5-lb. 
Figure  5-1  e  shows  the  subtracted  image  obtained  at  the  output 
plane  of  the  proposed  processor.  Since  we  have  utilized  an 
incoherent  source,  the  coherent  artifact  noise  in  the  subtracted 
image  is  substantially  suppressed. 

Another  interesting  experimental  demonstration  for  inspec- 
1 8 

tion  of  IC  mask  utilizing  a  simple  white- light  processor  is 
also  shown  in  Figure  5-2.  Since  the  white-light  source  con¬ 
tains  all  the  visible  wavelengths,  it  is  possible  to  identify 
the  defects  or  faulty  cracks  through  color  spatial  filtering. 

To  demonstrate  the  defect  detection,  we  assume  that  the 
Fourier  spectra  distribution  of  the  IC  mask  is  known  a  priori  . 

We  note  that  the  scale  of  t'ae  signal  spectrum  is  proportioned 
to  the  wavelength  of  the  light  source,  a  red  and  green  color 


Figure  5-2.  A  Wti  1 1  e -1. 1  g  h  t  Optical  Processor  for  IC 
Mask  Inspection. 
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spatial  filter,  as  shown  in  Figure  5-3a  is  utilized  in  our 
axpriment.  Since  the  faulty  cracks  are  generally  high  spatial 
frequency  signals,  they  would  be  expected  to  be  diffracted 
outside  the  red  color  region  of  Figure  5-3a,  Thus  the  faulty 
cracks  would  form  green  color  images  superimposing  with  red 
color  IC  mask,  as  shown  in  Figure  5-3b.  From  this  figure, 
we  see  that  the  faulty  cracks  can  be  easily  identified  through 
the  color  coded  image. 

In  concluding  this  section,  we  would  like  to  point  out 
that,  it  is  feasible  to  utilize  an  incoherent  and  white-light 
optical  processing  technqiue  for  micro-circuit  inspection. 

The  techniques  are  capable  of  improving  the  reliability, 
efficiency  and  cost  for  large  scale  automatic  precision  in¬ 
spection.  They  are  particularly  suitable  for  mi c r o -c i r c u i t 
board  inspection  and  robotic  assembly. 

V I .  Color  Hologram  Generation 
6 . 1  Introduction 

A  number  of  methods  have  been  investigated  for  storing 

color  images  holographically  on  black-and-white  film.  The 

color  image  is  then  retrieved  by  either  viewing  the  hologram 

19-21 

at  the  proper  angle  when  illuminated  by  white-light 
or  by  placing  it  in  a  white-light  optical  processor  where 

2  2-23 

the  appropriate  color  and  spatial  filtering  are  performed”  "  . 

3  e  c  a  u  s  e  the  light  needed  to  form  the  hologram  is  coherent, 
the  hologram  and  thus  the  reconstructed  image  are  corrupted 
with  coherent  artifact  noise.  If  the  holographic  fringes  were 
formed  with  a  point-pair  coherence  between  the  object  and 


a  s 
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reference  beams,  derived  from  extended  incoherent  sources,  the 

coherent  noise  would  be  suppressed.  In  this  section,  we 

explore  the  use  of  source  encoding  to  establish  this  point- 

3 

pair  coherence  requirement  and  demonstrate  a  system  that 
uses  it  to  form  a  hologram  of  a  2-D  color  transparency.  The 
color  image  is  retrieved  by  placing  the  hologram  in  a 
white-light  optical  processor  where  spatial  filtering  is 
p  e  r  f  o  r  ned  . 


o . 2  Color  Hologram  Construction 

A  system  for  constructing  a  hologra  hie  image  of  a  color 
transparency  is  shown  in  Fig.  6-1.  For  simplicity,  only  two 
primary  colors  (.  i .  e  .  .  ted  and  green)  are  used,  for  the  demon¬ 
stration.  This  system  is  similar  to  the  s ou r c e- enc o d ing  image 
subtraction  system  as  discussed  in  section  1.2  except  the 
gratings  are  displaced  at  a  distance  i  from  the  Fourier  plane 

is  a  color 


At  t  h  e  input  plane  P ,  ,  the  object  0,.  , 

1  J  1  ( x ,  y ) 


2  (  x  ,  ) 


the 
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be  bene  parallel  to  the  optical  axis  with  the  centers  dis- 

h 

placed  a  distance  of  h.  =  +  — —  i  respectively.  These  displace- 

1  _  kh  t 

ler.  .s  introduce  two  phase  factors  of  ?  =  +  = —  from  the 

1  f  *■ 

object  and  reference  beams  respectively.  The  corresponding 

complex  amplitude  distribution  around  the  optical  axis  at 

the  output  plane  would  be, 

_ ■  khl  .  kh l 

g(x,y)  =  Oj^Cx.y)  e  1  f2  x  +  0Q(x,y)  e1  f2  x.  (6-1) 

If  we  assumed  that  0  ^  (  x  ,  y  )  and  0  (x,y)  are  real  functions,  then 

a  hologram  would  be  formed  with  spatial  carrier  frequency  of 
2  kh  l 

P  =  - ^ —  in  the  x  direction. 

°  f2 

However,  for  the  two  beams  to  form  interference  fringes, 
a  mutual  coherence  requirement  must  exist  between  the  points 
of  the  two  beams  that  overlap.  To  establish  this  coherence 
requirement,  a  source  encoding  technique  is  used  which  estab¬ 
lishes  a  coherence  function  at  the  input  plane  of 
sin  ( Nt ' m4) 

f(x)  =  - sine  (  *t  ,  (6-2) 

N  sin(n  H)  Xt 

\  t 

where  s  is  the  width  of  the  slits  of  the  source  encoding 
mask,  N  is  the  number  of  the  encoding  slits,  d  is  the  spacing 


between  the  slits,  and  \  is  the  wavelength  of  the  extended 

Light  source.  Thus  by  choosing  d  =  ,  we  shall  obtain  the 

o 

proper  point-pair  coherence  requirement  for  forming  the  holo¬ 


graphic  fringes.  We  shall  also  note  that  the  spatial  frequency 


of  the  phase  grating  should  be  identical  to  the  spatial 


frequency  of  the  encoding  mask.  Thus,  separate  encoding  masks 
and  phase  gratings  are  used  for  each  primary  color  of  light 
sources  as  shown  in  Fig.  6-1. 


64 


[ 


Note:  if  the  two  encoding  masks  are  separated  by  a 

distanc.  equal  to  the  width  of  the  transparencies  O^(x.y) 
and  0  (x,y),  the  red  and  green  light  beams  would  be  separ¬ 
ated  at  plane  P  ^  and  would  be  passed  through  each  of  the 
gratings  respectively.  The  achromatic  transform  lens 

will  then  shear  the  two  images  0,  (x  y)  and  0  (x,y)  back 

i  o 

together,  such  that  the  red  and  green  images  will  be 
superimposed  around  the  optical  axis  at  the  output  plane  P, 
to  form  a  spatial  multiplexed  color  hologram. 

6 . 3  Reconstruction  of  Color  Image 

By  placing  the  recorded  hologram  in  the  white-light 
optical  processor  shown  in  Fig.  6-2,  a  color  image  would  be 
obtained  by  simple  spatial  filtering.  This  process  for 
reconstruction  of  the  color  hologram  image  uses  the 
dispersion  of  light  into  different  color  Fourier  spectra 
in  the  fourier  plane.  For  a  given  spatial  frequency  com¬ 
ponent,  the  displacement  of  the  smeared  fourier  spectra  in 
the  fourier  plane  is  proportional  to  the  wavelength  of  the 
light  and  the  spatial  frequency  of  the  input  object. 

The  spatial  frequencies  recorded  on  the  hologram  are 
inversely  proportional  to  the  wavelength  of  light  which 
passed  through  that  area  in  the  original  transparency 
when  the  hologram  was  constructed.  Thus,  for  example:  red 

light  passing  through  an  area  of  the  hologram  that 
corresponds  to  red  on  the  original  transparency  would  be 
diffracted  to  the  same  place  in  the  Fourier  plane  of  the 
white-light  optical  processor  as  green  light  passing  through 


an  area,  on  the  hologram  that  corresponds  to  green  information 


t  i  ca 
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on  the  original  color  transparency.  By  placing  a  narrow  slit 
over  this  region  of  the  smeared  Fourier  spectrum  a  properly 
colored  image  of  the  original  transparency  may  be  viewed  at 
the  output  plane. 

6 . 4  Experimental  Result 

In  our  experimental  setup  a  mercury  arc  lamp  with  a  green 

O 

interference  filter  centered  at  wavelength  5460A  was  used  as 
the  green  source.  A  HeNe  laser  which  was  passed  through  spinning 
ground  glass  was  used  as  the  red  color  light  source.  The  cor¬ 
responding  encoding  masks  have  dimensions;  for  s  =  2.9  pm, 
d  -  29  pm,  and  M  =  100,  for  the  red  encoding  mask  and  for 
s  =  2.5  pm,  and  N  =  100,  for  the  green  encoding  mask.  The 
spatial  period  of  the  phase  gratings  were  29  um  and  25  u-m  for 
the  red  and  green  gratings  respectively.  The  gratings  were 
placed  at  approximately  l  =  f,  the  focal  length  of  the  transform 
lens.  Thus  a  spatial  frequency  of  twice  that  of  the  corresponding 
phase  grating  would.be  encoded.  For  example,  spatial  frequencies 
of  34.5  i/mm  for  the  red  and  40  1/ mm  for  green  were  recorded. 
Figure  3  shows  the  input  object  0  ^  (  x  ,  y  )  used  to  construct  the 
color  hologram.  Figure  6-4  is  the  result  obtained  for  the 
reconstructed  color  hologram  image.  We  note  that  by  the  noise 
present  in  the  color  hologram  image  may  be  reduced  by  placing 
the  input  objects  in  a  liquid  gate. 

1 . 5  Remark 

By  using  an  adaptation  of  the  source  encoded  image  sub¬ 
traction  system,  we  were  able  to  construct  a  2-dimensional 
color  hologram.  Such  a  hologram  offers  easy  white-light 
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Figure  6-3. 


A  Black-and-White  Photograph  of  Original 
Color  Transparency. 


Figure  6-4. 


A  Black-and-White  photograph  of  Recon¬ 
structed  Color  Image. 
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reconstruction  of  the  color  image  and  the  possibility  of 
reducing  the  coherent  artifact  noise  during  the  construction 
process.  Although  the  system  was  described  and  demonstrated 
with  2  primary  colors  it  could  also  be  extended  to  three 
colors.  We  would  also  note  that  the  technique  may  also  be 
possible  to  generate  a  3-dimensional  holographic  image  for 
phase-type  object  of  definite  depth. 

VII.  Conclusion 

This  study  has  led  to  certain  definite  conclusions. 

In  particular  : 

1.  Incoherent  image  subtraction  processors  offer 

a  higher  signal-to-noise  ratio  as  compared  with 
a  coherent  image  subtraction  processor. 

2.  The  incoherent  image  subtraction  processor  is 
suitable  for  the  application  of  tracking  and 
identi fiation  for  moving  target. 

3.  It  is  feasible  of  utilizing  a  white-light 
source  for  image  subtraction  operation,  a 
technique  is  currently  processing. 

4.  The  white-light  image  subtraction  processor 

is  capable  of  performing  color  image  subtraction. 

5.  The  image  subtraction  technique  can  also  be 
utilized  for  precision  industrial  inspection, 
which  would  benefit  the  U.S.  Army  inspection 
technology  and  her  supporting  industries. 

6.  We  have  demonstrated  that  2-D  color  holograms 
can  be  generated  by  image  subtraction  processors 


with  extended  incoherent  sources. 
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7.  The  white-light  image  subtraction  technique  is 
capable  of  performing  in  a  real-time  mode,  a 
research  program  is  currently  investigating. 

The  net  effect  of  these  conclusions  is  to  emphasize 
the  truth  of  our  initial  assumption  of  the  white-light  image 
subtraction  processor  can  be  an  essential  part  of  the 
missile  tracking  and  identification  application.  This 
research  effort  in  real-time  tracking  program  should  be 
encouraged  to  pursuing. 
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